Universal non-linear conductivity near to an itinerant-electron ferromagnetic 

quantum critical point 
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We study the conductivity in itinerant-electron systems near to a magnetic quantum critical 
point. We show that, for a class of geometries, the universal power-law dependence of resistivity 
upon temperature may be reflected in a universal non-linear conductivity; when a strong electric 
field is applied, the resulting current has a universal power-law dependence upon the applied electric 
field. For a system with thermal equilibrium current proportional to T" and dynamical exponent 

z, we find a non-linear resistivity proportional to _E^(i+°)-i . 
PACS numbers: 72.10.Di, 72.90.+y, 75.30.Kz, 75.40.-s 



The notion of quantum criticality provides one of the 
few unifying theoretical principles of strongly correlated 
electrons [l[ 0; S]- It describes a range of phenomena 
in systems that are near to continuous, zero-temperature 
phase transitions; phase transitions that are driven by 
quantum rather than thermal fluctuations. In thermal 
equilibrium, quantum critical systems show character- 
istic spatial and temporal scaling in their response to 
external probes. For example, the conductivity of an 
itinerant-electron system near to a magnetic quantum 
phase transition has a power-law dependence upon tem- 
perature j5., %i . 

The behaviour of quantum critical systems out of ther- 
mal equilibrium has begun to attract growing attention 
over the past few years. Near to a quantum phase tran- 
sition all of the intrinsic energy scales of a system, other 
than the Fermi energy, renormalize to zero. In thermal 
equilibrium, the only remaining energy scale is the tem- 
perature itself. Because of this, quantum critical systems 
are particularly susceptible to being driven out of equilib- 
rium by external probes. In certain situations the univer- 
sal temporal scaling near to the quantum critical point 
may reveal itself in universal features of the steady-state 
adopted out of equilibrium; the out-of-equilibrium state 
being largely determined by a system's dynamics. 

Several recent works have addressed the question of 
whether universality persists when a quantum critical 
system is driven out of thermal equilibrium by the a; 
plication of a strong electric field. In particular, Refs 
and [i, 0] considered two-dimensional superconductor- 
insulator transitions , ^V\\ , the former in the case where 
the quantum dynamics and phase transition were con- 
trolled by coupling to a Caldeira-Leggett bath and the 
latter in the case of intrinsic superconducting dynam- 
ics. These systems can indeed display universality out of 
equilibrium 12 1 in both their current response [3, Q and 
their current noise statistics @. The triumph of Refs. 
and was to provide field-theoretical derivations of the 
scaling predicted by naive dimensional analysis. Numer- 
ical studies of related one-dimensional systems produced 



similar results [13|. 

Whilst these works provide interesting proofs of prin- 
ciple — and indeed, may yet be compared with experi- 
ment — most quantum critical systems that are studied 
experimentally are of a rather different type. The criti- 
cal modes at the superconductor-insulator transition are 
charged and couple directly to the electric field. A more 
typical situation has critical modes without a charge — 
often magnetic — which affect transport by scattering 
from electrons. Here we address the question of whether 
universal non-linear response in transport occurs in this 
more general setting. 

We find that given certain conditions on size and geom- 
etry, quantum critical itinerant magnets show a universal 
non-linear current response. For a long, narrow sample, 
with an electric field applied along its length, we predict 
a universal non-linear scaling of current with field given 
by 
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(1) 



where the thermal equilibrium resistivity is proportional 
to T" and z is the dynamical exponent. In the case 
of the Moriya-Hertz-Millis [1, 0, 0| model of the critical 
ferromagnet, a = (d + z — l)/z. Provided that certain 
constraints upon the dimension of the system are satis- 
fied, this result does not depend further upon the system 
dimensions. In the following, we will take some time to 
discuss this matter and compare our results to those of 
related works. 

We hope that these results will provide an alternative 
experimental window upon quantum criticality. Despite 
its successes, the theory of itinerant-electron quantum 
criticality has some puzzling problems; although power 
law dependencies upon temperature are seen experimen- 
tally, there is often a discrepancy between the observed 
and predicted powers in transport. Non-linear response 
may help to resolve this issue by providing two con- 
sistency checks: whether the equilibrium exponents are 
consistent with the out-of-equilibrium exponents through 
Eq.([T]); and whether the out-of-equilibrium exponent is 
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consistent with the Moriya-Hertz-MiUis theory. 

Our paper is outhned as foUows: we begin in Section 
1 with a general description of our scheme, paying par- 
ticular attention to matters of geometry and heat flows 
within the system. This will enable a heuristic derivation 
of our main results and a detailed comparison with the 
complementary work of Mitra et al. [l5| . In Section 2, we 
will begin with a survey of the Boltzmann treatment of 
the linear response of the itinerant-electron quantum crit- 
ical system in thermal equilibrium. This will allow us to 
introduce some notation and familiarize the reader with 
its application in this context. We follow this in Section 
3 by applying the Boltzmann transport formalism to the 
out of equilibrium system. This section contains a formal 
derivation of our main results. Finally, in Section 4, we 
turn to a discussion of the limitations of our analysis and 
of the prospects for seeing the effects that we predict in 
experiment. 

GENERAL SCHEME 

Geometry: We consider a long, narrow, quantum-critical, 
itinerant electron system with an electric field applied 
along its length. The system is longer than its transport 
length so that an electron traversing the sample scatters 
from magnons many times. The system must also be 
wide enough that it displays bulk behaviour, but narrow 
enough that heat generated within the sample can be 
transported to the boundary. 

Critical Fields: Starting from some low base tempera- 
ture. To, and gradually increasing the electric field one 
may anticipate two fields at which the response may be- 
come non-linear: 

i. When the energy gained by an electron from the elec- 
tric field between scattering events exceeds the tempera- 
ture; 



where Itr is the transport scattering length, 
ii. When the Joule heating rate exceeds the rate at which 
heat may be transported from the sample by a transverse 
heat flow; 



where a and k are the electrical and thermal conductiv- 
ities. The latter result has been obtained using k/ctTo = 
hh/hr- hh is the thermal scattering length and W is the 
sample width. 

In this work, we will be primarily concerned with the 
former case. In order for a system to be in this regime, 
we require that Ei ^ E2, so that we hit the field Ei 
first when increasing the electric field from zero; i.e. we 
require that W <ti \/hrhh- In addition, for the system 



to exhibit bulk behaviour requires that it be wider than 
its correlation length, W ^ Ith- Combining these two 
conditions upon the sample width yields 

Zth < < VU^- (2) 

Due to additional angular factors, the transport length is 
substantially greater than the thermal scattering length, 
^tr ^ ^th, so that there is a large window of sample widths 
over which the type of non-linear response that we en- 
visage can occur. In the high-temperature limit (with 
T <C ep nevertheless) in which experimental investiga- 
tions of itinerant electron quantum criticality are usually 
carried out, Ztr ~ Ith/S"^ where 9 ^ q/kp ~ (T/ep)^^^ is 
the angle of scattering. 

This regime is somewhat delicately balanced between 
the macro- and microscopic. In a truly macroscopic sam- 
ple where — > 00, non-linearity always occurs due to the 
failure to conduct away excess Joule heat. In our case, 
the transverse size of the system must be small enough 
that W <ti \/hihh, but the system inherits its behaviour 
from macroscopic equilibrium properties since it is larger 
than the correlation length 14]. If the above constraints 
are satisfied, the non-linear transport properties depend 
only upon bulk properties and are independent of the 
dimensions of the system. 

Thermal Coupling: Determining the non-linear response 
requires keeping careful track of the various heat fiows. 
We consider a simplified scheme of thermal couplings in 
our sample: 

i. The electrons couple to a heat sink at the boundaries 
of the sample and scatter from magnons. We do not con- 
sider electron-electron scattering since this is higher or- 
der in temperature or electric field than electron-magnon 
scattering and so sub-leading at low temperatures and 
fields. 

ii. The magnons may scatter both from one another and 
from the electrons. We do not consider coupling between 
magnons and the heat sink. Our reason is that magnon- 
phonon relaxation is higher order in temperature or field 
than magnon-electron scattering and therefore weaker at 
low temperatures and field. 

Heuristic Treatment: Given these descriptions of the ge- 
ometry of our system and the various microscopic cou- 
plings, we are now in a position to give a heuristic deriva- 
tion of our main results. Heat enters the system via Joule 
heating and ultimately leaves through a transverse heat 
fiow maintained by a transverse variation in tempera- 
ture. In the absence of scattering between electrons, this 
energy must pass through the magnons: Joule heating 
pumps energy into higher moments of the electron dis- 
tribution. This is ultimately carried away by a trans- 
verse heat flow maintained by a gradient in the symmet- 
rical part of the electron distribution. Energy can only 
pass into the symmetrical part of the distribution due to 
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scattering via magnons. In leading approximation, the 
magnons are raised to an effective temperature Tcff{E) 
determined by a balance between the Joule heating rate 
and the rate at which the magnons at Tcs{E) lose energy 
to the electrons at Tq ~ 0. Equating these two rates of 
change of energy leads to a self-consistency equation for 
Toff. 

In the high-temperature limit, the scattering time t, 
the transport scattering time rtr and the magnon energy 
decay rate d£^/dt are related as follows: 



1 

Ttr 

d£ 



r 

r2 



Using these relations for a system whose thermal equilib- 
rium resistivity scales as T", we find 

Joule Heating cx aE'^ oc E'^T'^ , 



Energy Relaxation oc t'^^ i/z+a)^ 
dt 



(3) 



By equating these two rates leads we deduce that Tcff oc 
^z/(z(i+a)-i) g^jjj j ^ ^j^^ following 

sections will flesh out these ideas with particular refer- 
ence to the Moriya-Hertz-Millis 0, M, m model of the 
critical ferromagnet. 

Comparison with Mitra et al: A recent work of Mitra et 
al has considered the same system as studied here, but 
in a different geometry. The results obtained in Ref.fisj 
are different from ours because of this geometry. In order 
to allay any confusion, it is worth spending a moment to 
note the main distinction between our two works. Mitra 
et a^[l5| consider an itinerant electron system with es- 
sentially two-dimensional geometry and an electric field 
applied in the third short direction. In this case, an elec- 
tron traversing the sample from one lead to another does 
not scatter appreciably from magnons — the electron dis- 
tribution is determined to leading order by the distribu- 
tions in the leads and may be written directly in terms of 
them using a Keldysh formalism. Mitra et al present an 
appealing derivation of this zeroth order distribution and 
show, using a renormalization group analysis, that an ef- 
fective temperature proportional to the applied voltage 
results. In our case, by contrast, an electron traversing 
the sample between the two leads scatters many times off 
magnons and the electron distribution must be calculated 
self-consistently from the start [l6l|. 

In the rest of this paper, we will spend some time flesh- 
ing out the mathematical details of this general scheme. 
We begin in the next section by reviewing the Boltzmann 
approach to thermal equilibrium transport in quantum 
critical metals. 



BOLTZMANN APPROACH IN THERMAL 
EQUILIBRIUM 

We will use Boltzmann transport techniques to analyse 
the out-of-equilibrium response of a quantum critical sys- 
tem to an electric field. Although this approach is famil- 
iar in other contexts, itinerant-electron quantum-critical 
transport is usually analysed by other means. Therefore, 
in this section, we will spend a little time summarising 
quantum critical transport in thermal equilibrium and 
how this may be described using a Boltzmann equation 
approach. This exposition will also serve as a useful way 
of defining the notation that we will use later in our anal- 
ysis of the non-linear response. 

Our first step will be to describe the thermal- 
equilibrium magnon propagator. We follow this by writ- 
ing down the electron Boltzmann equation and construct 
its linear response solution. Finally, we quote a number 
of relaxation rates that will useful in our non-equilibrium 
analysis. The details of the calculation of these within 
a Boltzmann framework is somewhat similar to that of 
the relaxation rates due to phonon scattering. We sketch 
these calculations in Appendix A. 



Magnon Propagator 

We work within the Moriya-Hertz-Millis [1, H, @] ap- 
proach to itinerant electron quantum criticality. The 
bosonic, magnetic critical modes — the magnons — are 
treated separately from the electrons (although they are, 
of course, made from electrons). The effects of scatter- 
ing between the magnons and electrons are treated self- 
consistently; the magnon dynamics being determined by 
Landau damping and the electronic transport being de- 
termined by scattering from the magnon fluctuations. 

The first step in the Moriya-Hertz-Millis approach to 
itinerant electron quantum criticality, is to determine 
the critical properties of the magnons. These critical 
properties are the combined result of the magnons' self- 
interaction and their overdamped dynamics due to Lan- 
dau damping. The simplest way to do this is through the 
self-consistent renormalization group Alternatively, 
one may use a more rigorous application of the renormal- 
ization group in order to obtain essentially the same 
results. In either case, the critical magnon propagator 
takes on the following form in the equilibrium quantum 
critical state: 



i^^ + 9' + r{T) 



(4) 



where Fg describes the Landau damping. Tq is propor- 
tional to |q| in the ferromagnet and constant in an anti- 
ferromagnet (or Fg = F|q|^~^ in general). The magnon 
gap r(T) takes on characteristic power-law forms in tem- 
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perature in the quantum critical system; 

r(T) oc T^^, 



(5) 



where d is the dimension and z is the dynamical exponent 
= 3 in the ferromagnet and 2 in the antiferromagnet). 
Through most of our subsequent analysis, we shall con- 
centrate upon the situation in three dimensions. This is 
readily extended to other dimensions. The overdamping 
of magnons has important consequences. Unlike their 
phonon counterparts, magnon excitations do not have a 
well-defined energy for a particular wave-vector. This 
does not have enormous consequences for a Boltzmann 
analysis in thermal equilibrium, but it does necessitate 
modification of the magnon Boltzmann equation when 
we consider the out-of-equilibrium situation. 



The Boltzmann Equation 

The electronic Boltzmann equation for scattering from 
an uncharged auxiliary mode may be written 



[dt + (eE/h) ■ dk] fk 
dq 



{2nf 



Ikqfk (1 - fq) - Iqkfq (1 " fk)] ■ 



(6) 



The matrices ^kq describe scattering from the auxiliary 
modes — magnons in our case. Quite generally, for scat- 
tering from auxiliary modes that have a thermal distribu- 
tion, the scattering matrices satisfy the detailed balance 
relationship 



Ikq = 7qfccxp[(efc - eq)/T] 



(7) 



In the case of magnon scattering, the scattering matrices 
take the form 



Ipq ~ \9q~p\ 



[1 + n{ep - eg)] p{p q,ep- e,) 
+n{eg - ep)p{q - p, Eg - ep) 



(8) 

where 5q_p is the matrix element for electron-magnon 
scattering and p{q,uj) is the magnon spectral function. 
In antiferromagnets, the matrix element Qq-p has signifi- 
cant momentum dependence, with scattering hot spots 
corresponding to resonance of the magnetic ordering 
wave- vector with the Fermi surface. For simplicity, we 
restrict our analysis to the case of ferromagnets or long 
wavelength helimagnets where the momentum depen- 
dence of 5q-p is weak and can be neglected. The magnon 
spectral function is given by 



1 



P{q.^) ^ --1rnD^ {q,uj) = — 

TT [r + q^Y + [uj/TqY 



(9) 



It is determined by the magnon propagator given in 
Eq. (ID). It contains all of the information about how 
dynamics is incorporated into the critical behaviour 



through the relative scaling of frequency and momen- 
tum: Lo ^ q'^ ^ ^^Tq. In what follows, it will prove very 
useful to work with the general form of the Boltzmann 
equation ^ rather than the form obtained after explicit 
substitution of 7pq. 



Linear Response Solution of the Boltzmann 
Equation 

The generic notation of Eq. ^ allows us to construct 
a formal linear response solution of the Boltzmann equa- 
tion both in thermal equilibrium and, ultimately, out of 
thermal equilibrium. In order to orient ourselves for the 
latter more involved calculation, let us first construct the 
conventional linear response solution with this general 
notation. Identifying 9] 



M.q = ^i^ Ik 

Ik I- fq 



dqikq^-^ (10) 

I - Jfc 



and adopting an Einstein convention with implied inte- 
gration over the momentum q, but not fc, we may write 
the Boltzmann equation in the form 

[dt + (eE/h) ■ dk] fk = -7fe[l - M]fcq/q(l - /q). (11) 

Let us consider an initial thermal distribution of electrons 
and auxiliary modes at the same temperature. The devi- 
ation in the electron distribution from its initial thermal 
distribution, fj^, in response to an electric field is given 
by a solution of the linearised equation 

{eE/h) ■ dk [f^ + 5fk] - -7fe[l - M]kqSfq, (12) 

where an Einstein convention has again been adopted. 
There are a couple of steps required in deriving this equa- 
tion. Firstly, we have used the fact that the scatter- 
ing integral is zero when the electrons and magnons are 
in thermal distributions at the same temperature. One 
must also allow for the dependence of [1— MJ^q upon fq in 
obtaining the first functional derivative of the scattering 
integral. 

A formal solution to the linearised Boltzmann equation 
(1121) is readily obtained. Expanding to linear order in the 
electrical field we find 



Sfk^[l-M]kq-E-dqfq. 



(13) 



This result may be integrated to obtain the current that 
flows in response to the application of the electric field; 



3 = 



dk 



dk dq 

(27r)3 (27r)3 



k[l-M]-^l-E-dqfq, (14) 



where an explicit integral over k has been restored. In 
the next section, we will turn to a consideration of the 
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non-linear response of the electron-magnon system using 
a very similar Boltzmann transport analysis. Before this, 
we identify a number of different time-scales of relevance 
to our problem and calculate them in the case of magnon 
scattering. 



A Compendium of Relaxation Rates 

The electronic scattering integral is given by the right- 
hand side of Eq. ^ or Eq. (fT2l) . 



dt 



dp 



lqpfq{^ - fp) 

-Ipqfpi^ - fq) 



scatt J c^^y 

= -lq[l - M]qpfp{l ^ fp 



(15) 



where jq and Mgp are given by Eq. ((10)) . Integration 
over p has been suppressed in the final two expressions, 
which are drawn from Eqs. pT|) and (HH), respectively. 

We may identify several different time-scales from this 
scattering integral that will appear in our later study of 
the non-equilibrium response; 



1 



1 

Q 

d£ 
'dt 



= lq = 



dp 
(2^ 
dp 



Iqp 



Iqp 



l-Jp 
fq 

-fp 



1 



1 



(27r)3 

dp dq 

(27r)3 (27r)3 



1-/. 



1 - 



g2 7*1- 



7qp/q"(l - /°) 

-7pq/°(l - /q°) 



(16) 



The first of these scattering rates is simply the inverse 
time between collisions. The second is the transport scat- 
tering rate. This has the usual additional geometrical 
factor arising since large angle scattering has more ef- 
fect upon transport than small angle scattering [l3|. The 
ratio 7q /7p is conventionally set to 1 since we are in- 
terested in the scattering of fermions near to the Fermi 
surface. The final expression is the rate of flow of energy 
from auxiliary modes at a temperature T (at which 7pq is 
evaluated) to electrons at temperature Tq (indicated by 
the superscript on the electron distribution functions). 

In the high-temperature limit, these relaxation rates 
have the following temperature dependence within the 
Moriya-Hertz-Millis [4i, |5|, |6,] theory in d dimensions: 

l/r cx T^d+z~3)/z^ 

d£/dt cx r(rf+3^-3)/^. 



(17) 



Details of how to get these results from Eqs. (|T6l) are given 
in Appendix A. After these preliminaries, we are now in 
a position to adapt our Boltzmann equation to describe 
the out-of-equilibrium behaviour of our system. 



NON-EQUILIBRIUM RESPONSE 

In this section we will turn the machinery of Boltzmann 
transport to the question of non-equilibrium behaviour in 
the itinerant critical ferromagnet. As discussed earlier, 
for a system to have an out-of-equilibrium steady state 
under the application of an electric field, it must be cou- 
pled to a heat sink that can dissipate the excess energy 
generated by Joule heating. We must pay careful at- 
tention to the various thermal couplings. The nature of 
these bears repetition at this juncture. 

As described in the introduction, we consider a long, 
narrow sample in which the excess Joule heat is carried 
away by a transverse heat current to a heat sink at the 
edge. Heat entering the electrons via Joule heating passes 
to the magnons and then back to the electrons via mu- 
tual scattering and leaves the electrons via coupling to 
a heat sink at the boundary. The magnons themselves 
interact both with the electrons and with one another. 
Electron-electron scattering is neglected in our analysis — 
it is higher order in temperature and hence field — as is 
coupling of magnons directly to the heat sink. 

Our analysis is divided into three parts. We be- 
gin by writing down the Boltzmann equations for the 
electron-magnon system. These equations embody the 
various thermal couplings and interactions in our sys- 
tem. The only subtlety enters through the form of the 
magnons' Boltzmann equation: the overdamped nature 
of the magnon excitations leads to a slightly more com- 
plicated equation than the comparable case of phonon 
scattering. In fact, the details of the magnon Boltzmann 
equation will not have a huge effect upon our main result. 
Next, we will present formal solutions for the electron and 
magnon distribution functions. Our main results follow 
from consideration of these solutions in the limit where 
magnon-magnon scattering leads to a thermal distribu- 
tion of magnons with temperature determined by the 
electric field. We will end with an argument why cor- 
rection to this thermal distribution of magnons do not 
change the scaling of our results. 



Tlie Boltzmann Equation 

The electron Boltzmann equation is given by a mini- 
mal modification of the thermal equilibrium Boltzmann 
equation ([6]): 



[dt + {eE/h) -dk+Vk-V] fk 
[f, n] + Scattering to heat sink 



I em 
'k 



dq 

(2^ 



Ikqfk (1 - fq) 

-Iqkfq (1 - fk) 



heat sink 



(18) 



where \'j^[f,n] indicates the scattering integral for elec- 
trons of momentum k scattering from magnons. The gra- 
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dient term w^.V has been added to allow for the possibil- 
ity of transverse heat flow. The electron-magnon scatter- 
ing matrices now take a slightly modified form compared 
with that in thermal equilibrium ([5]). Since the magnons 
are overdamped and as a result do not have a definite 
relationship between their energy and momentum, the 
magnon distribution is a function of both energy and 
momentum. Taking this into account, the scattering ma- 
trices take the form 



7pq l.9q— pi 



+nq-p{eq - ep)p{q -p,(q- ep) 



(19) 

The linearised expansion about the zero-field, base- 
temperature distribution /g takes the form 

{eE/h) ■ dk [fl + 5h] 



dq 



Ikqfk (1 - /q) 

-Iqkfq (1 - fk) 



(27r)3 

-7fe[l - M]kq5fq + heat sink 



(20) 



where 7fc and M^g take a slightly modified form out of 
equilibrium given by 



Ik 

7fcMfcq 



(|^[7fcq(l-/q")+7q;./q°] 



lkqfu+lqk{l-fu) 



(21) 



These reduce to our previous expressions for 7^ and M^g 
in thermal equilibrium (jlO[) . The simplified forms given 
by (fTO|) can be obtained by making use of the detailed bal- 
ance condition, which is not satisfied out of equilibrium. 
A couple of points are worth making about Eq. ipO)) . 
Firstly, there is a zeroth order term on the right-hand 
side. This term is not present in thermal equilibrium 
(it is zero upon applying the detailed balance condition). 
This term has a different symmetry in momentum space 
than the first order term in 5f and we will use this in our 
analysis shortly. 



The added complication due to the magnons being 
overdamped is compounded when we come to write down 
the magnon Boltzmann equation in a moment. Luckily, 
this does not affect the bulk of our calculation. We will 
use the formal notation ll™!/, n] through as much of our 
analysis as possible in order to keep algebra to a mini- 
mum. When we eventually substitute the particular form 
of the scattering integrals near to the end of the calcu- 
lation, we will find that most of the integration of these 
scattering integrals carries over directly from the thermal 
equilibrium calculation. 



The Magnon Boltzmann Equation takes the form 



dp dq 

(2^)3 (27r)3 '■^'^ 



~nk{e)fp{l - fq) 

-[l + nfe(e)](l-/p)/g 



S{e + ep - eq)6{p - q + k) 



J *i^^2A3^^^^^^p(ei,Pi)p(e2,P2)p(e3,P3) 



-nfc(e)np,(ei)[l-Knp,(e2)][l + np3(e3)] <:(i,,„ „ „ ^ ^ ^\ 
+ [1 + nk{e)][l + npjei)]np^{e2)np,{e3) \ 



(22) 



The easiest way to see the origins of the various terms 
in this equation is to momentarily treat the magnons as 
if they had a definite relationship between energy and 
momentum. In this case, the magnon spectral function 
becomes a delta-function and the scattering integrals re- 
duce to the same form as those for electron-phonon scat- 
tering. As for the electron Boltzmann equation, we will 
carry out as much of our analysis as possible using the 
formal expressions and for the magnon- 

electron and magnon-magnon scattering integrals. 



Solving the Boltzmann Equations 

Our analysis of the Boltzmann equations ([18]) and p2|) 
derived above proceeds as follows: We begin by divid- 
ing the electron distribution function into two parts — a 
spherically symmetric part and a non-symmetric part. 
The magnon Boltzmann equation is divided similarly. 
After this division, the resulting Boltzmann equations 
have simple interpretations. The equation for the sym- 
metric part of the distribution function describes the bal- 
ance between the transverse heat flow and the flow of 
energy out of the magnons into the symmetric part of 
the electron distribution. The equation for the remain- 
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ing part describes a balance between the flow of energy 
from the electrons into the magnons and the Joule heat- 
ing rate. In order to obtain useful results from these 
equations, we will go to a limit where the magnon distri- 
bution is assumed to be thermalised at some temperature 
TesiE) . The final step in our analysis will be to show that 
corrections to the thermal distribution of magnons do not 
change the way in which the current response scales with 
field. 

Expanding the Boltzmann Equation: 
The electron distribution function is divided into its sym- 
metric part Z*^ (assumed to be a Fermi distribution at a 
low base temperature that varies slowly across the sam- 
ple) and the remainder 5f. With this notation and after 
expanding to linear order in 5/, the Boltzmann equations 
may be written in the form 

v^-^fi = (23) 

J I em 

E-d^{fl + 5Q = + (24) 

^|me 

(25) 

We have adopted an Einstein convention where terms like 



{S\^ /5fk)5fk are implicitly integrated over k. The su- 
perscripts S and A refer to symmetric and non-symmetric 
parts of the scattering integrals in q. We have allowed for 
a transverse gradient in which supports a transverse 
heat flow. 

One might question how, given the fact that we are 
interested in the non-linear response, we can use a lin- 
ear analysis in 6f. In a linear-response, relaxation-time 
approximation, the Fermi surface is effectively shifted a 
distance Tn-eE/h in momentum space. Provided this is 
much less than the Fermi wavevector {rt^eE/h <^ fcp) a- 
linear response analysis may be applied. In the present 
case, it turns out that the transport relaxation time, rtr, 
self-consistently becomes a power of E so that the resul- 
tant current is non-linear in E. 

Heat flows 

The physical content of the equations ((23l) , ((24| and ((25)) 
is most readily appreciated by considering the energy 
transfers that they represent. In the case of ([23]) and 
(|24p we multiply by the electron energy Cq and integrate 
over q. In the case of (|25p . we multiply by the magnon 
energy e and the spectral density p{k, e) and integrate 
over k and e. After doing this, Eqs. (|23ll25l) reduce to 



= 
= 
= 



dq 
dq 
dk 



dep{e, k)e 



ime r j'O 
'k,eU 



Sfg 



(26) 
(27) 
(28) 



Ea. ([26)l may be interpreted as a balance between the 
transverse heat flow — described by the first term on the 
right-hand side — and the energy flowing into the sym- 
metrical part of the electron distribution — described by 
the second term on the right hand side. The flow of heat 
into the heat sink has been treated as a boundary condi- 
tion in writing down this equation. Solving this equation 
leads to the explicit limit on the sample width discussed 
in the introduction and worked out in detail in Ref. 
We will not concentrate upon it further here. 

Eg. ([27)1 can be interpreted as a balance between Joule 
heating — described by the first term on the right-hand- 
side — and the rate at which energy flows from the non- 
symmetrical part of the electron distribution into the 
magnons — described by the second and third terms. 
To see this requires a little manipulation. The first term 



may be written explicitly as Joule heating after integrat- 
ing by parts with respect to q. 

Ea. (p8)) corresponds to a balance between the rate at 
which energy flows into the magnons from the symmet- 
rical and non-symmetrical parts of the electron distribu- 
tion — described by the first and second terms respec- 
tively. The net flow of energy into the magnons from the 
electrons is zero in a steady state, if we neglect heat flow 
directly from the magnons to the heat sink. The latter 
process is ignored since it is much slower than magnon- 
electron scattering. The third term is identically zero, 
since magnon-magnon scattering conserves energy. This 
fact is extremely useful. By considering this integrated 
equation, one can avoid having to deal explicitly with the 
magnon-magnon scattering integral. 

We can transform this flnal equation into a more useful 
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form by using the fact that electron-magnon scattering is energy conserving. This imphes that 



J 



dq 5r^"'[f,n] 



(27r) 



3 eg 



Sfk - 



dk 



(27r) 

r 



gde ep(e,fc) — '-r? ^fq 



i.e the energy entering the electrons from the magnons 
is equal to the energy entering the magnons from the 
electrons. Using these results reduces Eq. ((28|) to the form 



dq 

(2^ 



(Sl^™[/°,nl 

ir[/^^]+ \i 

^ OJk 



(29) 



To make further progress, we must solve the Boltzmann 
equations explicitly. We will do this to leading order 
through an approximation that the magnon distribution 
is thermal. We will then argue that corrections to this 
do not alter the scaling. 

Thermal Magnon Approximation 

Our leading approximation is to assume a thermal distri- 
bution of magnons, 7ifc(e) ~ n'^{e) — (e*^/^""' — 1)^^, where 
the effective temperature is to be determined shortly. In 
this case, the linearised Boltzmann equations 
and (|^ reduce to 



N (5l""[/o,nO] 



5fk 



(30) 
(31) 
(32) 



We have used the fact that the magnon-magnon scat- 
tering is identically zero for a thermal distribution of 
magnons. The second equation may be formally solved 
for Sf to obtain, 



Sf 



E.da 



5f 



(33) 

where we have suppressed momentum labels and inte- 
grals over momentum for brevity. Expressions such as 
((5l'^™/(5/)~^ are to be understood as matrix inverses with 
appropriate integrations over momentum in their prod- 
ucts. In order to determine the effective temperature, we 
substitute this solution for Sf into the energy-integrated 
form of Eq. ((^ or (PTjl given by Eq. and expand to 



leading order in E. The result of this substitution is 

-1 







dq 



3 eg 



I em 



E-da 



6\' 



Sf 



E ■ d„Sf 



(34) 

Integrating the second term by parts reduces it to the 
form 



= 



dq 



3^1 



1 9 



SI" 



Sf 



dgSf 



(35) 



The second term is now explicitly the leading order con- 
tribution to the Joule heating rate. This is balanced 
against the first term which describes the decay of en- 
ergy from a thermal distribution of magnons at tempera- 
ture Tcff{E). Since both the electron and magnon distri- 
butions involved in the expressions are thermal distribu- 
tions — at T = and T = Tcs{E) respectively — we may 
evaluate Ea. ([35| using the results of Section 3. Writing 
Eq. ([35|) in terms of the scattering matrices of Section 3, 
it can be reduced to 



dq 



(27r)3 

dq 
(2^ 



eg7g[l-M]gp/p(r0)(l-/p(r0)), 

3eg7g[l-Af]gp<5/p. (36) 



Substituting for 5f to leading order in E from Eq. ([ST 
into Eq. ([5^ . we obtain 



dq 

(2^ 



■eg.ag[[i-M],pi7pi£;.ag/g(ro)] 



dq 

(2^ 



Joule Heating 
eq7q[l-M]gp/p(To)[l-/p(ro)] (37) 



Energy decay from T^ato To 

This equation may be written in the form d£/dt cx E^Ttr 
as before in Eq,([3]). Since the magnon distribution is 
thermal at temperature Toff, using the temperature scal- 
ing of the various relaxation rates given in Ea. (|17p . we 
find 
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in the high-field/temperature limit and d dimensions, im- 
plying a non-linear current 

j cx £:(^-i)/('i+2z-2) 



as suggested in Section 1. 

Corrections to Thermal Magnon Approximation 

Corrections to a thermal distribution of magnons may 
be rather large, since in the absence of the electric field 
the magnons are essentially in a zero-temperature distri- 
bution. We argue, nevertheless, that corrections to the 
thermal distribution of magnons considered above do not 
change the scaling of current response. The analysis is 
similar to the calculation of phonon drag in thermal equi- 
librium (which similarly does not change the scaling with 
temperature) . 

We expand the magnon distribution to linear order 
about the effective thermal distribution; n = nP -\- 6n. 



Substituting into the Boltzmann Eqs. (l23|24p and (|25l) . 
we obtain 



E ■ d, (/o + Sfg) 



ir[/^»°]+ \ (38) 



Sfk 



-5fk 



Sn 



Snt, 



k.e 



= 



(39) 

"^'[y ' Sfg (40) 



5n 



Sn 



^q,«- 



In Eas. (|38)) and ([39]), 6n has been divided into symmetric 
and non-symmetric parts 6n^ and 6n^. These contribute 
to the equations for the symmetric and non-symmetric 
parts of the electron distribution respectively. Eq. lHTI) 
can be solved formally for Sn with the result 



Sn 



^|me ^|mm 



5n 6n 



5n^ 



in 



Sf 



(42) 



(43) 



We have identified the spherically symmetric and non- 
symmetric parts of Sn. Substituting this back into 



Eq.(|38D and ^ one obtains 



§\em /rime rimmX "1 



Sn \ Sn 
S\T^[f,n°' 



Sfk 



5fk 



Snk,e 



Sn 



sr 



SV 



Sf 



■Sf 



(45) 

In the second of these equations, we have adopted the no- 
tation of Lifshitz-Pitaevskii^l^ identifying this as a term 
describing a magnon-mediated electron-electron interac- 
tion. The argument that magnon drag does not affect 
scaling is completed by showing that S\'^"^'^/Sf scales with 
at least as high a power of T as S\™^ /Sf. This requires us 
to go beyond the generic form of the scattering integrals 
to use their explicit expressions for magnon-electron scat- 
tering given in the Boltzmann equations (fT8|) and (|22p . 
We ignore the magnon-magnon scattering (it is higher 
order in T — and hence E — than the magnon-electron 
scattering) we may write 



5I« 



Sf 



Snk, 



SV 



SV 



Sn 



SV 



Sf 



(46) 



Taking the explicit form of the scattering integrals, the 
various functional derivatives may be written as 



Snk^e 

s\r. 



\9?p{k,e) 



ifq - fk-q)S{(- - Cq + Efc-g) 



+ {fq - fk+q)S{e - Eq+fc + Eq) 



(41) Sni^, 



= \g\H{l - k)S{v - e) 



dp 

(2^ 



(/p+fc - fp)S{<i + - Ep+fe) 



Sfi 



\9\ 



-(?^fec + fi+k)S{e + ei- ei+k) 

+ (1 + nkc - fi-k)S{e + ei-k - q) 



Substituting these equations back into Ea. (|46)) shows 
that the corrections due to magnon drag lead to con- 
tributions to the electron scattering integral that are at 
least of the same order in temperature as the direct con- 
tribution. 



CONCLUSIONS AND PROSPECTS 

We have considered non-linear transport near to an 
itinerant electron quantum critical point. Since the dy- 
namics near to a quantum critical point are universal, 
and since steady-state, out-of-equilibrium distributions 
are determined by dynamics, we have argued that the 
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universality present near to an equilibrium quantum crit- 
ical point may be reflected in the out-of-equilibrium be- 
haviour. 

There are two ways in which a quantum critical itiner- 
ant electron system may be driven out of equilibrium by 
an electric field. At the highest fields, the rate of Joule 
heating overwhelms the rate at which heat may be trans- 
ported out of the system by thermal conduction and the 
system heats up until the two balance. This mechanism 
leads to non- linear response in truly bulk samples. We 
have considered non-linear response in a restricted ge- 
ometry where we anticipate that conductivity becomes 
non-linear at a lower field governed by the rate at which 
energy can scatter between electrons and magnons. The 
resulting conductivity is expected to be independent of 
sample size and geometry (provided that certain con- 
straints are satisfied). At the lowest fields, the response 
will return to the linear, thermal equilibrium response. 

The existence of the intermediate range of non- 
linearity requires restrictions upon the sample width so 
that thermal conduction can be maintained at a sufficient 
rate to transport away heat generated by Joule heating. 
The sample width must nevertheless be sufficient that the 
magnons demonstrate their bulk behaviour- i.e it must 
be larger than the magnon correlation length. Because 
of the rather different scaling of transport and thermal 
relaxation lengths with temperature (and hence field), 
there is a large window of fields within which the type of 
non-linearity that we investigate should exist. 

What are the prospects for seeing these effects exper- 
imentally? We have described a particular experimental 
geometry in which the heat current is transverse to the 
electrical current. This enabled the algebra to be readily 
negotiated. In order to see these effects experimentally, 
we suggest a slightly different geometry [22|. One pos- 
sibility is the following: take a bow tie shaped sample 
with current injected and removed along opposite wings 
of the bow tie. A current sent through this sample should 
demonstrate a non-linear steady-state of the type that 
we have described. The constriction at the centre of the 
bow tie will have enhanced field and current densities 
and will operate in a non-linear regime. The injection 
of current along the extended edge of the bow tie will 
reduce contact heating; performing the experiment in a 
pulsed manner will further mitigate these effects. The 
large heat capacity of the wings of the bow tie compared 
to the constriction will allow a relatively long pulse time 
before heat effects become significant; the wings will ef- 
fectively act as a low-temperature heat sink. A sketch of 
this arrangement is shown in Fig. 1. 

It remains to estimate what field strengths give rise 
to the non-linear effects that we anticipate. For a typical 
quantum-critical itinerant magnet (e.g. Sr3Ru207, which 
has n = 2 X 10^^ m^'^, cr = lO^il^^m^^) and typical 
cryogenic temperatures of around 100 mK, the electric 
field required to observe these non-linear effects is of the 



Non-linear Response region 




FIG. 1: Schematic diagram of proposed experimental 
system: i. Current enters and leaves the bow tie shaped 
sample along the edges of the wings, reducing contact heating, 
ii. Enhanced field and current density in the constriction leads 
to non-hnear response in this regime, iii. The extended wings 
act as a low-temperature heat sink. 

order 

£;=M:«50Vm-i 
evpT 

(We have used the Drude formula a = ne^r/m and n = 
47rfcp/3) or a voltage drop of 0.25 V over a typical sample 
length of 5 mm. 

In conclusion, the universal power-law scaling of con- 
ductivity near to an itinerant, magnetic, quantum criti- 
cal point is reflected in a universal power law scaling with 
electric field in the non-linear conductivity regime. This 
provides a new way to investigate the consistency be- 
tween theoretically predicted power-laws and those seen 
experimentally. 
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APPENDIX A 

In this appendix, we outline the evaluation of the scat- 
tering rates given in Eq. (fT6|) . The details of the calcula- 
tion is somewhat similar to that of the same relaxation 
rate due to phonon scattering. Because of this similar- 
ity, our analysis follows quite closely — and makes ex- 
plicit reference to — the calculation of electron-phonon 
relaxation rates presented in Chapter 8 of the textbook 
of Mahanji^. We shall carry out the calculations in d- 
dimcnsions, where d = 2 or 3 in the physical system. 

Our first task is to make a few manipulations of the 
energy relaxation rate to put it in a simplified form. The 
first step involves using the detailed balance relation ex- 



11 



pressed in the form 



Ipq — Iqp rrp 



/J) 



- /J) 



where n^{Aeq 
and Aeo = Ca 



is the Bose distribution at temperature T 
- Cp. The next step is to integrate over |qr| 
assuming that the |g|-dependence of terms other than 
the Fermi-distribution functions is small and also that 
the density of electronic states is constant at the Fermi 
surface. In so doing, we encounter two integrals 



/iH = j def{e){l- f{e-uj)) = ujn{uj), 
hH = / dej{e - - /(e)) = c^(n"(cc>) + 1). 



After these manipulations, the energy relaxation may be 
expressed as 

where pp is the electronic density of states at the Fermi 

surface and the remaining integral over q is just angular; 

functions of q are to be interpreted as having |qr| equal 

to the Fermi wavevector. The superscripts and T on 

the Bose-distribution function indicate that they are at 

the base temperature or the elevated temperature, T, 

respectively. Notice that this is automatically zero when 
^0 _ 2^ 

In the limit T° 0, we may neglect n'^(Aeq). Also, 
in the limit where Ae^ ^ T, we may write the energy 
relaxation in the same form as the other relaxation rates 
using n^(Ae,)/(n^(A6q) + l) ^ (1-/J)/(1-/J). As aU 
the terms in our scattering rates are now at the elevated 
temperature T, we will from now on omit this superscript 
for clarity. 

So far, we have reduced our scattering rates to the 
forms 



l-/p 





f d'p 








f d'^p 


T-tr ' Q 

q 




d£ 1 / 


d'^p d'^q 


H ~ ~2 J 


{2-kY (27r)'i 



l-/q 

l-/p 
l-/q 



1 - 



gp7q 

q2 



PFAe^7qp 



l-/p 
l-/q 



(47) 



In the present case, we are interested in scattering from 
critical magnons. 7pq then takes the form given by 
Eq. ([5]) . In order to calculate the various relaxation rates 
explicitly, it is useful to introduce a generalization of the 
McMillan function. In the discussion of electron-phonon 
scattering, this takes the form(23j 



(48) 

where tOq is the frequency of a phonon with momentum 
q and Cp is the energy of an electron with momentum p. 
The generalisation of this to scattering from overdamped 
modes is given by 



a'FiE,u;)^^ J ^|.g,pp(q, ^)^(£; _ efc+,), (49) 



where p{q,uj) — T'mD^{q,uj) is the magnon spectral 
function. The analogous McMillan function for trans- 
port is given by 



(50) 

The additional angular factors weight scattering at dif- 
ferent angles in the usual way. Eqs. (|49l50p are the gen- 
eralisations of Eqs.(8. 145,8.146) in Ref.^]. 

With this identification, expressions for the various 
scattering integrals may be obtained in the limit Aeq <C 
T as follows (see for example Section 8.3.1 of Ref.[23| — 
we use (1 - /p)/(l - fq) ~ n(Aeq)): 
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1 27r f°° f°° cfiq , 

— = 2— / dujn{u})a^F{tk,uj) ^2 I div——^\gqfn{Lj)p{q,uj)6{ek+q-ek), (51) 

Tfe I Jo Jo l^Tl'J 

1 2tt f°° f°° d'^q q k 

= 2— j dujn{ui)a'^F{ek,uj) = -2 J dw-^^--p|gqpn(w)-^p(q', w)5(efe+q - efe), (52) 

d£ f°° f°° d'^q 

— = AttH^pf / dujn{uj)uj'^a'^F{ek,uj) = 2h^pF / duj \gq\^n{Lj)uj'^ p{q, uj)5{€k+q - £fc). (53) 

Jo Jo (.^""j 



These integrals may be simplified by first linearising the 
electron energy near to the Fermi surface; ek+q ^ Cfc ~ 
Vk Q = vp cos9\q\, where Vk is the Fermi velocity and 9 
is the angle between k and q. Assuming that the matrix 
element does not have a significant angular dependence — 
an assumption that is only true for ferromagnets; anti- 
ferromagnets have hot lines of scattering on the Fermi 
surface where electron states are related by the order- 
ing wave-vector that lead to complications in this latter 



case 
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the angular integrals over q may then be car- 
ried out. One then obtains 



1 

Tfe 



2g' 



(27r)2«F Jo 



dujn(uj) 



X / d\q\\q\^-'piq,u;) 



2g^ 1 
(27r)2«F 2fcf Jo 



dujn{uj) 



X j d\q\\qfpiq,u;) 



d£ 



2g 



2 r°° 



dujn{ijj)Lo'^ 



{2t:?v^ Jo 

T/vp 

X / dlqWql'-^piq^iu) 



(54) 



The momentum integral has acquired an explicit cut-off 
at T/v-p after linearising the electron energy at the Fermi 
surface. The remaining integrals may be calculated after 
explicit substitution of the magnon spectral function. 

We will evaluate these expressions in both high- and 
low-temperature limits. Which of these limits one is in 
is determined by a comparison of r{T) with the typical 
value of The former varies with temperature accord- 
ing to r{T) ^ and the latter as T^. In the low- 
temperature limit, r ^ q^ and in the high-temperature 
limit r <^ q^. In both cases, we consider temperatures 
much less than the Fermi energy, T <C ep. 

Carrying out the integrals in the low-temperature limit. 



we find 



1 



rd^nico) d\q\\q\^~^p{q,u) 
[2'kYvy Jo Jo 



d\q\\qf'^ 
Jo 



duj nluj}— -. — 



d\q\\qr'Tq —n{LU 

Jo '^J^ q 



duj 



Uj/Vq 



"I rT„'^'~'l + {u/rTqY 



d\q\\qV-^Tq dun{urTq)- 



T 

W)Jo 



d\q\\q 





d-2 



du 



1 + U2 



where we have used the fact that r ^ q^ a.t low tem- 
peratures. The frequency integral is dominated by the 
region where uj is up to order r(T)rq. For q ^ T, at 
these frequencies w/T ~ Tr{T)T^~^ <C 1 and the Bose 
distribution can be approximated by its low-frequency 
limit n{x) ~ T/x. As a final consistency check, we need 
to make sure that the dominant momentum q is of the 
order of T as indeed it is. 



A similar evaluation of the transport scattering rate 
yields the result 



T 



'd+2 



r{T) 



(55) 
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The energy relaxation is given by 



given by 



2„2 ^oo 

h^PF / dujuj^n{uj) I 'd\q\\q\'^~'^p{q,uj) 



{2^?VF Jo 



r2 + (u/Tg)^ 



tT r^qWqr^Tl 



T tt' 
^ 2 



/ "'d\q\\q\''-^T^ 



TM d\q\\q 
Jo 

rpd+Z—l 



d+z- 



(56) 



d£ 



25^ 



X / ' d\q\\q\''-'p{q,ui) 



rf-"(-Vy^ ^l9lrw-^ + i 



dLon(uj)uj 



dv n{vT)v 



d+z-4 



du- 



u^^ + 1 

(59) 



Carrying out the same integrations in the high- 
temperature hmit, we find 



1 

Tfc 



rdc.n(c.) /"d|q||q|'^-V(q,c.) 
(27r)2z;F Jo Jo 



Lo/T\q 



z-2 



q^ + iu;/r\q\^-^y 



"'o y 



n{uj) 



|d+z-4 



(ia;n((jj)w 





d-3 



"'0 
/>oo 

d+s-3 / , ^ d-3 



U"" + 1 

_ i/Zi+^-i 

T ' dv n{vT)v ' du 

Jo Jo u 



2z 



1 



(57) 



In carrying out these manipulations we have used the fact [12] 
that ^ ^ oo at high temperatures, which is con- 

sistent since the dominant contribution to the frequency 
integral comes from a; ~ T. We have rescaled the mo- 
mentum and frequency integrals and then and used the 
explicit substitution Tq — V\qY~'^ . A similar evaluation 
of the transport scattering rate yields 



(58) 



i.e. it carries an extra factor of q^ ~ T"^/^ compared with 
the scattering rate. Finally, the energy relaxation rate is 
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